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Abstract 



The pathway idea is a way of going from one family of functions to another family of functions and yet 
another family of functions through a parameter in the model so that a switching mechanism is introduced 
into the model through a parameter. The advantage of the idea is that the model can cover the ideal or 
stable situation in a physical situation as well as cover the unstable neighborhoods or move from unstable 
neighborhoods to the stable situation. The basic idea is illustrated for the real scalar case here and its 
connections to topics in astrophysics and non-extensive statistical mechanics, namely superstatistics and 
Tsallis statistics, Mittag-Leffler models, hypergeometric functions and generalized special functions such as 
the H-function etc are pointed out. The pathway idea is available for the real and complex rectangular 
matrix variate cases but only the real scalar case is illustrated here. 

1. Introduction 

Mathematical techniques in the area of special functions, statistical techniques in the area of statistical 
distribution theory and characterizations and information theory techniques in the area of generalizations 
of Shannon type entropies and their axiomatic definitions and properties had been developed by the first 
author in the period of time from 1965 to 1980 [Mathai and Rathie (1975), Mathai and Pederzoli (1977), 
Mathai and Saxena (1978)]. Then starting from the 1980's applications of these techniques to astrophysical 
problems were explored by the present authors jointly in the areas of energy generation, solar modelling, 
and gravitational instability problems. A summary of the findings until 1988 is available from Mathai and 
Haubold (1988) and later results from Mathai (1993), Mathai and Haubold (2008), and Mathai, Saxena, and 
Haubold (2010). It was seen that interesting results, mathematically and statistically and with potential of 
physical interpretations, could be obtained by the fusion of special function theory, statistical distribution 
theory, and information theory. The present authors' work in basic space sciences started in the 1980's and 
then in 1988 a sequence of United Nations Workshops, prospectively co-organized by ESA, NASA and JAXA 
and hosted by individual nation states was conceptualized at the Centre for Mathematical Sciences, India 
(http:/ /www. cmsintl.org/). 

For new results in mathematics and mathematical statistics, applications are usually found in sooner 
or later interplay with natural sciences and hence the vast amount of theoretical development that has been 
achieved, do not all have immediate applications. One sequence of results which are useful in pathway 
model building in physical situations based on experimentally recorded data was demonstrated recently 
[Haubold et al. (2012)] . An idea was introduced in the 1970's by which one could go from one family of 
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functions of mathematical or statistical nature to another family to yet another family, and later in 2005 
[see Mathai(2005), Mathai and Haubold (2007)] the idea was extended to cover real and complex scalar 
mathematical or random variables, as well as real and complex rectangular matrix variables. The basic 
notion will be explained with the help of a specific example. Consider a general input-output type situation. 
It could be reactions generating or annihilating particles, diffusion and transport of particles and thereby 
an emerging product is what is observed. Consider particle reactions and let N(t) be the number density at 
time t and the rate of reaction denoted by [Mathai and Haubold (1988)]. If the number of particles 

produced is proportional to the original population size then the differential equation is d ^*- ) = A N(i) 
where A denotes the rate of reactions [Saxena et al. (2010)]. Let the diffusion rate or destruction rate be \x 
then the residual rate is c = A — \i. If production dominates then c > and if destruction dominates then 
c < 0. Then for the model 

- -c N(t) => N(t) = N e- ct (1.1) 

where iV is the initial population size. If the rate of change is proportional to a power of the population 
size and if decay dominates then the equation and the solution are the following: 

±N(t) = -c[N(t)F => N(t) = -[1 - c(l - a)t]T^. (1.2) 

This is a power law type of behavior. For a < 1 the function in (1.2) belongs to a particular case of a type-1 
beta family of functions. Let N(i) in (1.2) be denoted by N\(t). For a > 1, by writing 1 — a = —(a — 1) 
and denoting N(t) by N 2 (t), we have 

JV 2 (t) = [l + c(a-l)t]-^r. (1.3) 

Here (1.3) is a special case of a type-2 beta family of functions. When a — > 1, denoting N(t) by A^i) in 
this case, 

N 3 (t) = lim N 2 (t) = lim N^t) = c~ ct . (1.4) 

t — t — ^ 1 _ 

This, in fact, is the model in (1.1). Thus, N\(t) and N 2 (t) for a < 1 and a > 1 respectively describe a wide 
range of models. If the exponential form in (1.1) is the stable form in a physical situation then a here can 
be called the stability parameter and Ni(t) and N 2 (t) can describe the unstable neighborhoods of N 3 (t). 

2. Optimization of Entropy 

Models in physical situations are also constructed by optimizing entropy measures. The basic Shannon 
entropy in a probability scheme, for the continuous situation is 

/oo 
f(x)\nf(x)dx (2.1) 
-oo 

where f(x) is a statistical density and k is a constant. When k is present, we can assume f(x) to be any 
non-negative integrable function, need not be a statistical density. S represents a measure of uncertainty in 
a probability scheme. If S(t) is maximized over all functional / satisfying the condition f(x)dx = 1 
and f{x) > for all x then / is the uniform density. If (2.1) is maximized subject to two conditions (i): 
/-oo /( x ) a - x = 1 an d E(x) is a given quantity, E{x) = J_ x f(x)dx — the expected value or the mean 
value of x then we end up with / being an exponential density. With reference to the number density in 
(1.1)- (1-3) the second condition will imply that the expected number, E[N(t)] in a unity space in unit time 
is a fixed quantity which can also be interpreted as the principle of conservation of energy if we are dealing 
with energy generation. If, further, the second moment E(x 2 ) is also fixed then we end up with the Gaussian 
or normal density. 

The basic entropy measure in (2.1) is generalized in different directions. A class of a-generalized en- 
tropies, their characterizations and properties are available from the book: Mathai and Rathie (1975). One 
of the a-generalizcd entropies, in the continuous case is 

[f°° (f(x)) 2 - a dx- 1] 
M a (f)= °° Q -l (3-2) 
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Consider the optimization of (3.2) subject to the conditions 

/CO 
\x\ s f(x)dx = k\ < oo, 
-co 

/CO 
\x\' y+5 f(x)dx = k 2 < oo 
-oo 

where k\ and fc 2 are fixed, and the optimization is done over all non-negative integrable functions. 7 = 0,5 = 1 
is the case leading to (1.1) to (1.3) or Tsallis statistics [Tsallis (1988)]. Consider the function g(f) over all 
functional /, where 

9(f) - lf(x)?- a Ar \x\\f(x) + A 2 \xP +s f(x) 
where Ai and A 2 are Lagrangian multipliers. Then the Euler equation is given by 

^jg(f) = (2 - a)[f(x)] l - a A \x[< + A 2 \x\^ +s = 

=> /(z) = dl^Hl - o(l - a)\x\ 5 ]^ (3.3) 

where is taken as c\ and ^ i s taken as a(l — a), a > 0. Note that (3.3) for a < 1, a > 0, 5 > 0, x > 
can be called an extended generalized type-1 beta model. For a > 1, writing 1 — a = —(a — 1), (3.3) reduces 
to the following: 

/ 2 (z) =c 2 |x| 7 [l + a(a-l)|a;| l5 ]-^r,a> l,5>0,a>0. (3.4) 

Note that (3.4) can be called an extended generalized type-2 beta model. Denoting f(x) under a < 1 as 
fi(x) we have 

h(x)= lim h{x)= lim / 2 (^) - c 3 |xre- a l-l 5 (3.5) 

which can be called an extended generalized gamma model. This is the entropic pathway. If fi(x), / 2 (x) 
of (3.3)-(3.5) are taken as statistical densities then c\, c 2 , C3 can act as the normalizing constants, which are 
available by integrating out in (3. 3), (3. 4) and (3.5) respectively. 

Cl = [a(1 ' Q)] ^ r(2 ff + T 7 +I) ,a<l,a>0,5>0, 7 +1>0 (3.6) 
c 2 = n : pr,a > 1 (3.7 

2 r(2+±)r(^-2±i)' 
1 7 + 1 

a > 0, S > 0, 7 + 1 > 0, - > and 

a — 1 

T+l 

c 3 = ,a>0,<5> 0,7+1 >0. (3.8) 

21 V 8 > 

The model in (3.3) for a general a is the scalar version of the pathway model of Mathai (2005). This is the 
distributional pathway. Here a is called the pathway parameter. When a < 1 then the model describes the 
whole family of functions belonging to extended generalized type-1 beta family. When a > 1 then we move 
into the whole family of functions belonging to the extended generalized type-2 beta family. When a — > 1 
then both these families go into the family of extended generalized gamma family. We can also look into the 
transitions in the corresponding differential equations. This is the differential pathway. Thus we have the 
following pathways [(Mathai and Haubold (2007)]: 

Entropic pathways 

Distributional pathways 

Differential pathways. 
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Note that (3.3) for x > 0, 7 = 0, S = 1, a = 1 is Tsallis statistics of non-extensive statistical mechanics 
which works for all the cases of a < l,a > l,a —¥ 1. This particular case of (3.3) is also the model in 
(1.2). It is said that over 5000 articles are produced on Tsallis statistics so far since it was introduced in 
1988 [Tsallis (1988), see also Hamza (2005)]. It is also said that over 3000 people are working on this model 
giving various types of interpretations in various fields. 

Model (3.3) for a > 1,6 = l,a = l,x > is what is known in the literature as superstatistics [Beck 
and Cohen (2003), Beck (2006)]. Note that since superstatistics assumes the functional form in (3.4) for 
a > 1, from superstatistics one cannot get (3.3) for a < 1. In the family of pathway models, superstatistics 
is derived from the case a > 1 and a — > 1 whereas Tsallis statistics covers all cases a<l,a>l,a— >1 but 
the main restriction here is that 7 = or the factor x 7 is absent in Tsallis model. In superstatistics x 7 is 
present but it covers only the type-2 beta (a > 1) and gamma (a — > 1) families of functions and not type-1 
beta (a < 1) families of functions. 

In the actual applications, x could be time, energy, velocity etc and then the models can have different 
interpretations in different disciplines. 

4. Bayesian Procedure 

The model in (3.5) for a prefixed parameter a can be written as a conditional density of the type 

h{x\a) = \xPe- a W\a > 0,-oo < x < 00. (4.1) 

Suppose that the parameter a has a prior density given by 

9(a) = ^f^« e ~ le " ? ' a > 0, r/ > 0, e > (4.2) 
where e and 77 are known constants. Then the unconditional density of x is given by 



/ 

J a 
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7 + 1 



f 4 (x\a)g(a)da = - J > +1 , / a^+^e-^+N >da 

) Ja=0 

< 5 l-( 2 T i +^) 



2?f r(e)r(7+±) Ja=0 

N 7r(l+l + e ) r l 



2 r(2+i)^r( e ) l v 



x 



|7r(2+i + e )q 
2r(2+i)r( e ) 



T+l 



j^—ll + V^^. (4.3) 



This (4.3) for x > is the superstatistics. Note that for the convergence of the integral i + |x| 5 must 
remain positive. Hence superstatistics can only produce type-2 beta family of functions when considering 
gamma type conditional density for x\a and gamma type marginal density for a. When r/ is of the form 
b(a — 1), b > 0, a > 1 and 7 -ti + e = — ^ then we have the pathway model for a > 1. The unconditional 
density of x in (4.3), denoted by f x {x), can also be interpreted the following way: fi(x\a) is the density 
of x where a is a parameter. Then we are superimposing another density g(a) on the density fi(x\a) 
and then the resulting density f x (x) can be called superimposed statistics or superstatistics. Apparently 
when superstatistics was introduced they were unaware of Bayesian procedures in Probability/Statistics. In 
Bayesian procedure, superstatistics is the unconditional density of x when x and the parameter a, for which 
a prior density is assumed, both belong to gamma family of densities. A more general family of unconditional 
densities is available from Mathai and Haubold (2007). Dozens of papers are published on superstatistics 
and it is being hotly pursued in different disciplines. 

5. Fractional Considerations 

Going back to our basic growth-decay problem where the rate of change is proportional to the population 
size, our basic differential equation, equation (1.1), is 

j t f(t) = -c f{t), c> => f(t) - f = -cj f(t)dt. (5.1) 
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If the total integral is replaced by a fractional integral of the Ricmann-Liouville type let us see what happens. 
The left sided Riemann-Liouville fractional integral operator is denoted by oD~ a = qI" and it is defined as 

oD- a f = -^- ( X ( x -ty- 1 f(t)&tMu)>Q- (5.2) 

"(a) Jo 

Fractional integral can be given many interpretations in statistical literature as fraction of a total integral, as 
the density of residual variable u = x — y where x and y are independently distributed real positive random 
variables such that x — y > etc [Mathai (2010)]. If the total integral in (5.1) is replaced by fractional 
integral of (5.2) then the equation becomes 

f(x)-.f = -c( D- a f)(x) (5.3) 

where /o is a constant. One simple way of solving this equation is by taking Laplace transforms on both 
sides. Let the Laplace parameter be s. Let the Laplace transform of / be denoted by f(s). Then 

Lf(s) fo / e~'*dx = -c / e— [-— / (x - t)^ 1 f(t)dt]dx. 

Jo Jx=0 r («) Jo 



Then 



/ - T = - s " a /» / = n /" - a] (5-4) 
s s[l + cs a \ 



oo 

C 



= /oE(^) fe (- 1 ) fe - 

k=0 

Taking the inverse Laplace transform we have 

f(x) = h^i-lt^^ = foEai-af) (5.5) 



where E a (-) is the basic Mittag-Leffier function. Generalization of the basic Mittag-Leffler function are the 
following: 



E ^ x ) = E T(l + ak) Ma) > °' El{x) = ° X 



,_„-(! + «*) 

OO 

E <*>?W = E T (/+ak) Ma) > > ° 

KA') = E^fc8T^)' K(a) > °'^ (/3) > (5 ' 6) 

where (7)^ is the Pochhammer symbol 

(7)fc = 7(7 + l)-(7 + k-l), (7)0 = 1,7^0. 

More generalized form of (5.5) is the Wright's function, which is a special case of the H- function. More on 
the applications of these functions may be seen from Mathai and Haubold (2008), Mathai et al. (2010). 

It is seen that when we move from a total differential equation to a fractional differential equation, 
Mittag-Leffler function and its generalizations, Wright function and H-function enter into the solutions. A 
series of recent papers are available on the solutions of fractional reaction equations and fractional reaction- 
diffusion equations. The Laplace transform in (5.4) belongs to a general class of Laplace transforms, see 
Mathai et al. (2006) and the various references therein, and various members from this general class appear 
when solving some fractional differential equations. Some of the papers may be seen from Haubold et al. 
(2011) and Saxena et al (2010). 
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The effects of power transformations and exponentiation on various models can be seen from a recent 
paper Mathai (2012). Let us see what happens if a parameter is becoming larger and larger in a Mittag- 
Lcffler model of (5.6). Suppose that j3 is real and it is becoming larger and larger. Then by using the 
asymptotic expansion of gamma functions or as a first approximation the Stirling's formula 

T(z + a) ss V27rz z+a ~5 C ~ z for \z\ -» oo, a is bounded 

we see that 

a k (j) k x 5k v / 2^-i<=-' 3 



oo k f \ 

= E^«§) , )* = ( 1 + aa ' a )~ 7 - ( 5 - 7 ) 

fc=0 ' P 

This is the pathway model, Tsallis statistics and superstatistics for the case a > 1 for 7 = ^-j- , a = 
b(a — 1), b > 0, a > 1, 5 > 0. For 7 = 1, (5.7) becomes a power series. 

6. A Mathematical Perspective 

Mathematically speaking the whole process of transition from one functional form to another, Tsallis 
statistics, superstatistics and pathway models in the scalar case can be described as getting rid off some 
parameters from a hypergeometric series. Take for example a 1.F1 series: 

If we wish to get rid off an upper or lower parameter then we do a limiting process. 

x 

lim iFi(a;6; -) = ^i( \b]x) 

a— >oo a 

lim iF\{a\b]bx) = iF (a; ;x), \x\ < 1 

lim !Fo(a; ; -) = i J1 ( ; ;ar) = (6.2) 

a— 5-oc a 

lim Fi( ;b;bx) = F { ; ; a:) = e x . (6.3) 

The binomial function going to the exponential function in (6.2) is the basis for the pathway idea, Tsallis 
statistics and superstatistics. Observe that a similar rich class of pathways are available from (6.3) where a 
Bessel function is going to an exponential function. All the above limiting forms are available by using the 
fact that 

l im = i = i im ^L. (6.4) 

a— >oo a a— >oo [ajk 

All these ideas are extended to the matrix-variate cases, to real positive definite, hermitian positive definite 
and to rectangular matrices, see the basic paper Mathai (2005), and later papers by the author and his 
co-workers are also available. One such model is the following: 

f(X) = c \AiXBX'Aip\I-a(l-a)AiXBX'Ai\^ (6.4) 

where lisapxr, r >p matrix of full rank p of distinct real random or mathematical variables, A is a p x p 
constant positive definite matrix, B is a r x r constant positive definite matrix, X' denotes the transpose of 
X, A? denotes the positive definite square root of the positive definite matrix A, f(X) is a real-valued scalar 
function of X and c is a constant. This c can act as a normalizing constant if f(X) is treated as a statistical 
density. If the matrix X is relocated at some other matrix M then replace X by X — M in the model. The 
constants r\ > 0, a > and a are real scalars where a is the pathway parameter. For a < 1 the model in 
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(6.4) will stay in the generalized real matrix-variate type-1 beta family of functions. For a > 1 the model in 
(6.4) will go to the generalized real matrix-variate type-2 beta family of functions. When a — > 1 both these 
type-1 beta and type-2 beta families will go to a generalized matrix-variate gamma family of functions. This 
can be seen by using the result 

lim \I - a(l - a)A^XBX'Ai\^ = cxp{-an tr(A^ XBX'A^)} 

where tr(-) denotes the trace of (•). It can be seen that all the real matrix-variate densities that are used in 
the current literature are available from the model (6.4) for various values of the pathway parameter a. A 
similar rich family is there if we consider the transition from a Bessel form to the exponential form. Model, 
corresponding to the one in (6.4), is available when the variables are in the complex domain also. The results 
are parallel. 
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